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Abstract. The diffuse Galactic 7-ray emission is 
produced by cosmic rays (CRs) interacting with the 
interstellar gas and radiation field. Measurements 
by the Energetic Gamma-Ray Experiment Telescope 
(EGRET) instrument on the Compton Gamma-Ray 
Observatory (CGRO) indicated excess 7-ray emission 
> 1 GeV relative to diffuse Galactic 7-ray emis- 
sion models consistent with directly measured CR 
spectra (the so-called "EGRET GeV excess"). The 
excess emission was observed in all directions on 
the sky, and a variety of explanations have been 
proposed, including beyond-the-Standard-Model sce- 
narios like annihilating or decaying dark matter. 
The Large Area Telescope (LAT) instrument on the 
Fermi Gamma-ray Space Telescope has measured the 
diffuse 7-ray emission with unprecedented sensitivity 
and resolution. We report on LAT measurements 
of the diffuse 7-ray emission for energies 100 MeV 
to 10 GeV and Galactic latitudes 10° < |6| < 20°. 
The LAT spectrum for this region of the sky is well 
reproduced by the diffuse Galactic 7-ray emission 
models mentioned above and inconsistent with the 
EGRET GeV excess. 
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I. Introduction 

The diffuse 7-ray emission, both Galactic and ex- 
tragalactic, is of significant interest for astrophysics, 
particle physics, and cosmology. The diffuse Galactic 
emission (DGE) is produced by interactions of CRs, 
mainly protons and electrons, with the interstellar gas 
(via 7T° -production and Bremsstrahlung) and radiation 
field (via inverse Compton [IC] scattering). It is a direct 
probe of CR fluxes in distant locations, and may contain 
signatures of physics beyond the Standard Model, such 
as dark matter annihilation or decay. The DGE is a 
foreground for point-source detection and hence influ- 
ences the determination of their positions and fluxes. It 
is also a foreground for the much fainter extragalactic 
component, which is the sum of contributions from 
unresolved sources and truly diffuse emission, including 
any signatures of large scale structure formation, emis- 
sion produced by ultra-high-energy CRs interacting with 
relic photons, and many other processes (e.g., |Q] and 
references therein). Therefore, understanding the DGE 
is a necessary first step in all such studies. 

The excess diffuse emission > 1 GeV in the EGRET 
data Q relative to that expected from DGE models 



consistent with the directly measured CR nucleon and 
electron spectra [3| led to the proposal that this emission 
was the long-awaited signature of dark matter annihi- 
lation J4). More conventional interpretations included 
variations of CR spectra in the Galaxy 0, (3), contribu- 
tions by unresolved point sources [6|, and instrumental 
effects 0, H. 

A model of the DGE depends on the CR spectra 
throughout the Galaxy as well as the distribution of the 
target gas and interstellar radiation field (ISRF). Starting 
from the distribution of CR sources and particle injection 
spectra, the distribution of CRs throughout the Galaxy 
is determined taking into account relevant energy losses 
and gains, then the CR distributions are folded with the 
target distributions to calculate the DGE |9|. Defining 
the inputs and calculating the models are not trivial 
tasks and involve analysis of data from a broad range of 
astronomical and astroparticle instruments iflOll . 

The Fermi LAT was launched on June 11, 2008. It 
is over an order of magnitude more sensitive than its 
predecessor, EGRET, with a more stable response due 
to the lack of consumables. The LAT data permit more 
detailed studies of the DGE than have been possible ever 
before. 

In this paper, analysis and results for the DGE are 
shown for the Galactic mid-latitude range 10° < \b\ < 
20° measured by the LAT in the first 5 months of the 
science phase of the mission. This region was chosen 
for initial study since it maximises the fraction of signal 
from DGE produced within several kpc of the Sun 
and hence uncertainties associated with CR propagation, 
knowledge of the gas distribution, etc., should be min- 
imised. At lower Galactic latitudes the large-scale DGE 
is dominant while the emission at higher latitudes is 
more affected by contamination from charged particles 
misclassified as photons and uncertainties in the model 
used to estimate the DGE. Further details of the present 
analysis are given in ifTTI 

II. LAT Data Selection and Analysis 

The LAT is a pair-conversion telescope with a pre- 
cision tracker and calorimeter, each consisting of a 
4x4 array of 16 modules, a segmented anti-coincidence 
detector (ACD) that covers the tracker array, and a 
programmable trigger and data acquisition system. Full 
details of the instrument, onboard and ground data 
processing, and other mission-oriented support are given 
in flH. 
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The data selection used in this paper is made using 
the standard LAT ground processing and background 
rejection scheme. This consists of two basic parts: first a 
simple accept-or-reject selection (pre filtering) followed 
by a classification tree (CT) ifTSIl based determination of 
the relative probability of being background or signal. 
The prefiltering phase screens particles entering the LAT 
for their charge neutrality using the tracker and ACD. 
The direction reconstruction software extrapolates found 
particle trajectories in the tracker back to the scintillation 
tiles of the ACD covering the LAT and only accepts 
events in which the intersected tile shows no significant 
signal. In addition, considerations such as how well the 
found tracks project into the measured energy centroid 
in the calorimeter, and the shape of the shower energy 
deposition, are also used. The overall background rejec- 
tion of the prefiltering phase is 10 3 — 10 4 depending on 
energy, yielding a 7-ray efficiency > 90% relative to 
7-rays that convert in the LAT. 

Classification trees, which afford an efficient and 
statistically robust method for distinguishing signal from 
noise, are used to reduce backgrounds further. The CTs 
are trained using Monte Carlo data which have passed 
the prefiltering described above. Multiple CTs are used 
to make the procedure robust against statistical fluctu- 
ations during the training procedure. The result from 
averaging these CTs is the probability for an event to be 
a photon or background. This final rejection parameter 
allows the background levels to be set according to 
the needs of the analysis. For the analysis of diffuse 
emission, the CT generated probability is set to allow a 
Monte Carlo predicted orbit-averaged background rate of 
~ 0.1 Hz integrated over the full instrument acceptance 
> 100 MeV. This yields a 7-ray efficiency > 80%, 
and the remaining background is at a level where the 
majority of the contamination arises from irreducible 
sources such as 7-rays produced by CR interactions in 
the passive material outside the ACD, e.g., the thermal 
blanket and micrometeroid shield for the LAT. The 
events corresponding to the above criterion are termed 
"Diffuse" class and are the standard low-background 
event selection. 

The analysis presented here uses post-launch instru- 
ment response functions (IRFs). These take into account 
pile-up and accidental coincidence effects in the detector 
subsystems that are not considered in the definition of 
the pre-launch IRFs. Cosmic rays, primarily protons, 
pass through the LAT at a high rate and sufficiently 
near coincidences with 7-rays leave residual signals 
that can result in 7-rays being misclassified, particularly 
at energies < 300 MeV. The post-launch IRFs were 
derived using LAT events read out at regular intervals 
as a background overlay on the standard simulations 
of 7-rays and provide an accurate accounting for the 
instrumental pile-up and accidental coincidence effects. 
The on-axis effective area for the event selection used 
in this paper is ~ 7000 cm 2 at 1 GeV and is energy 
dependent; this is approximately 10% lower than the pre- 



launch effective area corresponding to the same event 
selection. The systematic uncertainties of the effective 
area, evaluated by comparing the efficiencies of analysis 
cuts for data and simulation of observations of Vela, are 
also energy dependent: 10% below 100 MeV, decreasing 
to 5% at 560 MeV, and increasing to 20% at 10 GeV 
and above. The point spread function (PSF) and energy 
resolution are as described in fl2l . 

The LAT nominally operates in a scanning mode that 
covers the whole sky every two orbits (i.e., 3 hrs). We 
use data taken in this mode from the commencement 
of scientific operations in mid-August 2008 to the end 
of December 2008. The data were prepared using the 
LAT Science Tools package, which is available from 
the Fermi Science Support Centre. Events satisfying the 
Diffuse class selection and coming from zenith angles 
< 105° (to greatly reduce the contribution by Earth 
albedo 7-rays) were used. To further reduce the effect 
of Earth albedo backgrounds, the time intervals when 
the Earth was appreciably within the field of view 
(specifically, when the centre of the field of view was 
more than 47° from the zenith) were excluded from this 
analysis. This leaves 9.83 Ms of total livetime in the 
data set. The energy-dependent exposure was calculated 
using the IRFs described above. 

The photon counts and exposure were further pro- 
cessed using the GaRDiAn package, part of a suite 
of tools we have developed to analyse the DGE fl4l . 
Gamma-ray skymaps with 5 bins per decade in energy 
from 100 MeV to 10 GeV were generated. For each 
energy bin the intensity was obtained by dividing the 
in-bin counts by the spectrally-weighted exposure over 
the bin. We used two methods for the spectral weighting: 
a power law with index —2 and the spectral shape of the 
assumed DGE model (described below). With the energy 
binning used in this paper the differences in the derived 
intensities were < 1% between these two weighting 
schemes. 

Figure Q] (left) shows the LAT data averaged over all 
Galactic longitudes and the latitude range 10° < |6| < 
20°. The hatched band surrounding the LAT data indi- 
cates the systematic uncertainty in the measurement due 
to the uncertainty in the effective area described above. 
Also shown are the EGRET data for the same region 
of sky derived from count maps and exposures available 
via the CGRO Science Support Centre [j and processed 
following the procedure described in O and we have 
included the standard systematic uncertainty of 13% 
fl5l . For both data sets the contribution by point sources 
has not been subtracted. The LAT-measured spectrum is 
significantly softer than the EGRET measurement with 
an integrated intensity Jlat(> 1 GeV) = 2.35 ±0.01 x 
10~ 6 cm~ 2 s _1 sr _1 compared to the EGRET integrated 
intensity J E gret(> 1 GeV) = 3.16 ± 0.05 x 10~ 6 
cm~ 2 s _1 sr _1 where the errors are statistical only. Not 
included in the figure is the systematic uncertainty in 

'http://heasarc.gsfc.nasa.gov/docs/cgro/egret/ 
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Fig. 1: Left: Preliminary diffuse emission intensity averaged over all Galactic longitudes for latitude range 10° < 
|&| < 20°. Data points: LAT, red dots; EGRET, blue crosses. Systematic uncertainties: LAT, red; EGRET, blue. 
Right: Preliminary LAT data with model, source, and UIB components for same sky region. Model (lines): ir°- 
decay, red; Bremsstrahlung, magenta; IC, green. Shaded/hatched regions: isotropic, grey/solid; source, blue/hatched; 
total (model + UIB + source), black/hatched. 



the energy scale, which is conservatively estimated from 
comparison between Monte Carlo and beam test data as 
< 5% for 100 MeV to 1 GeV, and < 7% above 1 GeV 
where it is believed that if any bias is present energies 
are overestimated. Taking the uncertainty on the energy 
scale into account, the LAT spectrum could be softer, 
increasing the discrepancy with the EGRET spectrum 
further. 

Figure Q] (right) compares the LAT spectrum with 
the spectra of an a priori DGE model, and a point- 
source contribution and unidentified background (UIB) 
component derived from fitting the LAT data that are 
described below. The DGE model is an updated version 
of the "conventional" model from GALPROP |9|. Major 
improvements include use of the formalism and corre- 
sponding code for pion production in pp-interactions by 
fl6l . a complete recalculation of the ISRF ifTTl . updated 
gas maps, and an improved line-of-sight integration 
routine. 

The source and UIB components were obtained by 
fitting the LAT data using GaRDiAn with the DGE 
model fixed. Point source locations were taken from the 
3 month Fermi LAT source list down to sources with 5- 
o significance. Due to the limited statistics of all but the 
very brightest sources, we used 3 bins per energy decade 
in the fitting procedure. Source positions were fixed but 
the spectra were fit using one free parameter per energy 
bin. The UIB component was determined by fitting the 
data and sources over all Galactic longitudes for the 
high-latitude region \b\ > 30° for the full LAT energy 
range shown in the figure. Using this high-latitude region 
minimises the effect of contamination by the bright 



Galactic ridge which can be significant even up to ~ 10° 
from the plane due to the long tails of the PSF at low 
energies. 

To determine the uncertainty of the source and UIB 
components, we modified the effective area to the ex- 
tremes of its systematic uncertainty defined before and 
refitted the data. Since the DGE model is fixed, the 
absolute change in intensity caused by the modification 
to the effective area propagates directly to the source 
and UIB components. The systematic uncertainty on 
these components is energy dependent and due to several 
effects. 

For energies > 10 GeV the PSF is - 0.2° (68% 
containment) and the sources are well-localised spatially. 
Since the model is fixed and the sky maps are sparser at 
high latitudes for the data taking period in this paper, the 
UIB component absorbs almost all of the intensity from 
the modification to the effective area. At low energies 
the PSF is wider, 3.5° (68% containment) at 100 MeV 
for 7-ray conversions in the front section of the LAT, and 
the sources are less well-localised spatially. In addition, 
the sky maps are well populated even at high latitudes 
and display spatial structure. The PSF broadening of the 
sources provides spatial structure and because the DGE 
model is fixed, more intensity is assigned to the source 
component to compensate in the fit. These effects lead 
to the systematic error in the source component being 
relatively larger than the isotropic at low energies and 
vice versa at high energies. Note, this applies for the 
high-latitude region from where the UIB component is 
derived, and also for the mid-latitude range for which 
we show the combined contribution by sources in Fig. Q] 
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(right). Because the uncertainties in the source and UIB 
components are not independent we have conservatively 
added their systematic uncertainties for the total intensity 
band shown in Fig. Q] (right). 

The UIB component comprises the true extragalac- 
tic diffuse 7-ray emission, emission from unresolved 
Galactic and extragalactic sources, and residual particle 
backgrounds (CRs that pass the 7-ray classification 
analysis and 7-rays produced by CR interactions in the 
passive material outside the ACD) in the LAT data. In 
addition, other relevant foreground components that are 
not completely modelled, such as emission from the 
solar disk and extended emission [18|, fl9l . ll20l and 
other potentially relevant "diffuse" sources |2TI . Il23ll . 
11221 are included. Hence, the UIB component does not 
constitute a measurement of the extragalactic diffuse 
emission. Furthermore, comparison with the EGRET 
estimate of the extragalactic diffuse emission ll24ll is 
problematic due to the different DGE models used and 
analysis details that are beyond the scope of the current 
paper and will be addressed in a subsequent publication. 

III. Discussion 

The intensity scales of the LAT and EGRET have 
been found to be different with the result that the 
LAT-measured spectra are softer. In our early study of 
the Vela spectrum 0251 . which was made using pre- 
launch IRFs, the difference was apparent already above 
1 GeV. Following on-orbit studies new IRFs have been 
developed to account for inefficiencies in the detection 
of 7-rays in the LAT due to pile-up and accidental 
coincidence effects in the detector subsystems. The 
inefficiency increases at lower energies, with the result 
that the IRFs used in the present analysis indicate greater 
intensities in the range below 1 GeV, with the magnitude 
of the effect ranging up to 30% at 100 MeV. Our 
confidence that the IRFs used in the present analysis 
accurately represent our knowledge of the instrument 
comes from detailed instrument simulations that were 
validated with beam tests of calibration units, and to 
post-launch refinements based on actual particle back- 
grounds. The systematic uncertainty on the effective area 
gives an energy dependent measure of our confidence in 
the IRFs used in the present analysis. 

As a consequence, the LAT-measured DGE spectrum 
averaged over all Galactic longitudes for the latitude 
range 10° < |6| < 20° is systematically softer than 
the EGRET-measured spectrum. The spectral shape is 
compatible with that of an a priori DGE model that 
is consistent with directly measured CR spectra. The 
excess emission above 1 GeV measured by EGRET is 
not seen by the LAT in this region of the sky. 

While the LAT spectral shape is consistent with the 
DGE model used in this paper, the overall model emis- 
sion is too low thus giving rise to a ~ 10 — 15% excess 
over the energy range 100 MeV to 10 GeV. However, the 
DGE model is based on pie-Fermi data and knowledge 
of the DGE. The difference between the model and data 



is of the same order as the uncertainty in the measured 
CR nuclei spectra at the relevant energies |26j. In ad- 
dition, other model parameters that can affect the 7-ray 
production rate (e.g., Xco- tne conversion between CO 
line intensity and molecular hydrogen column density 
in the interstellar medium) have not been modified in 
the present paper. Overall, the agreement between the 
LAT-measured spectrum and the model shows that the 
fundamental processes are consistent with our data, thus 
providing a solid basis for future work understanding the 
DGE. 
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